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ABSTRACT 



We present a new measurement of the scaling relation between X-ray luminos- 
ity and total mass for 17,000 galaxy clusters in the maxBCG cluster sample. Stack- 
ing sub-samples within fixed ranges of optical richness, A^ooj we measure the mean 
0.1-2.4 keV X-ray luminosity, (L x ), from the ROSAT All-Sky Survey. The mean 
mas s, (M 2n n), is measured from weak gravitational lensing of SDSS background galax- 
ies (| Johnston et al.l [20071 ). For 9 ^ A^oo < 200, the data are well fit by a power-law, 
(L Jf )/10 42 /i- 2 ergs- 1 = (12.6±\ A 3 (stat) ± 1.6 (sys)) ((M 200 }/10 14 h' 1 M ) Lffi±al3 . 
f- ^ The slope agrees to within 10% with previous estimates based on X-ray selected cata- 

logs, implying that the covariance in Lx and -/V200 at fixed halo mass is not large. The 
luminosit y intercept is 30%, or 2a, low er than determined from the X-ray flux-limited 
sample of iReiprich fc Bohringerl (|2002l ). assuming hydrostatic equilibrium. This slight 
difference could arise from a combination of Malmquist bias and/or systematic error 
' in hydrostatic mass estimates, both of which are expected. The intercept agrees with 

that derived bv lStanek et al.l (|2006[ ) using a model for the statistical correspondence 
between clusters and halos in a WMAP3 cosmology with power spectrum normaliza- 
tion as — 0.85. Similar exercises applied to future data sets will allow constraints on 
the covariance among optical and hot gas properties of clusters at fixed mass. 
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1 INTRODUCTION 



Bulk properties of galaxy clusters, such as galaxy richness 
and velocity dispersion, X-ray temperature and luminosity, 
mean weak gravitational lensing shear and X-ray hydrostatic 
mass, display strong internal correlations (for reviews, see 
iRosati et aiTl2002l : IVoitll2005l ). Although these correlations 
are anticipated by dimensional arguments — the "bigger 
things are bigger in all measures" perspective — detailed 
theoretical expectations are complicated by baryon evolu- 
tion uncertainties. Even with some baryon model prescribed, 
the stochastic nature of the halo evolution in a hierarchical 
clustering framework will impose variance about the mean 



behavior, the scale of which is likely to be sensitive to the 
baryon prescription. 

Improved understanding of cluster scaling relations is 
desirable on both astrophysical and cosmological grounds. 
Better astrophysical models of galaxy and hot gas evo- 
lution in halos will improve constraints on cosmological 
parameters obtained from studies of cluster samples (e .g. 
IVikhlinin et al.ll2006l : TRozo et al-lboOTl : iMantz et al.ll2007h . 

The relation between X-ray luminosity Lx and halo 
mass M is a diagnostic of the halo baryon fraction and the 
entropy structure of the intracluster gas. From X-ray data 
alone, constraints on the Lx-M relation have been obtained 
using either direct or statistical arguments. The direct ap- 
proach assumes that the intracluster gas is in hydrostatic 
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equilibrium, so that observations of the X-ray surface bright- 
ness, which provides Lx, and the X-ray temperature pro- 
file can be combined to estimate ther mal pressure gradients 
and, h ence, cluster binding masses. iReiprich fc Bohringerl 
(2002, hereafter RB02) measured the Lx~M relation in this 
manner for the X-ra y flux limited HIFLUG CS sample of 
roughly 100 clusters l|Finoguenov et al]|200ll ). The statisti- 
cal approach, on the other hand, assumes a form for the 
likelihood, P(Lx\M,z), which is convolved with the space 
density of halos in a given cosmology to predict cluster 
counts as a function o f X-ray flux and red shift. This ap- 
proach was applied by IStanek et aD (120061. hereafter S06) 
to the flux-limited REFLEX survey i|Bohringer et al.ll2004T ) 
using halo space densities calibrated for a ACDM cosmology 
l| Jenkins et al.ll200ll : lEvrard et al.ll2002T l. 

A third approach to the Lx~M relation employs weak 
gravitational lensing mass estimates derived from shear pat- 
terns of background g alaxies in deep optic al i maging of 
cluste r fields. Recently. iBardeau et al.l (|2007h and iHoekstral 
|2007l ) have performed such analysis of 11 and 20 clusters, 
respectively, at redshifts near z ~ 0.25. These studies pro- 
vide a powerful cross check between the two methods, but 
provide a weaker constraint on the mass scaling as they focus 
on the X-ray luminous tail of the cluster population. 

In this paper, we extend the dynamic range of the third 
approach using a large, uniform sample. We present the Lx~ 
M relation derived from ~ 17, 00 optically selected g alaxy 
clusters in the maxBCG catalog l|Koester et a l. 2007b). Bin- 
ning clusters by optical richness, A^oo, we measure the mean 



ROSAT All-Sky Survey (RASS: Voges 


Il999l). Details of 


the method are described in Rvkoff et al. 


(2007, henceforth 



R07). The mean mass, M200 = (-M200IA200IU of each bin 
is measured from analysis of the stacked weak gravitational 
lensing signal from SPSS galaxy images. Mea surements of 
the shear are p resented in ISheldon et all (|2007h . and masses 
are derived bv lJohnston et alT i 2007h . 

We begin with a brief description of the input data and 
the methods. We then present measurements of Lx and 
M200 for richness-binned maxBCG clusters, and compare 
the resultant scaling relation to the relations found by ap- 
plying direct or statistical arguments to X-ray flux-limited 
samples. Throughout the paper, we assume a flat, ACDM 
cosmology with Ho = 100 /ikms" 1 and Q, m = 1 — Qa = 0.3. 
In addition, our measurements are quoted at z = 0.25, close 
to the median redshift of the cluster catalog. Following R07, 
the X-ray luminosities are measured in the 0.1-2.4 keV ob- 
server frame, and fc-corrected to rest frame 0.1-2.4 keV at 
z = 0.25. 



within r2oo- Here, we briefly describe the maxBCG cluster 
catalog, binning strategies, and our methods for calculating 
mean mass and X-ray luminosity. 



2.1 maxBCG Catalog 

The maxBCG cluster catalog provides sky locations, red- 
shift estimates, and richness values for the cluster sample we 
employ. Details of the selection algorithm and catalog p rop- 
erties are published elsewhere l|Koester et al.l [2007al fbh. In 
brief, maxBCG selection relies on the observation that the 
galaxy population of rich clusters is dominated by bright, 
red galaxies clustered tightly in color (the E/S0 ridgeline). 
Since these galaxies have old, passively evolving stellar pop- 
ulations, their g — r color closely reflects their redshift. The 
brightest such red galaxy, typically located at the peak of 
galaxy density, defines the cluster center. 

The maxBCG catalog is approximately volume lim- 
ited in the redshift range 0.1 ^ z 0.3, with very ac- 
curate photometric redshifts (5z ~ 0.01). Studies of the 
maxBCG algorithm applied to mock SDSS catalogs indi- 
cate that the completenes s and purity are v ery high, above 
90% (|Koester et al.ll2007bl ; IRozo et a.l.ll2007l ). The maxBCG 
catalog has been used to investigate th e scaling of galaxy 
velocity dispersion with cluster richness (|Becker et al.ll2007r ) 
and to derive constraints on the powe r spectrum norma liza- 
tion, erg, from cluster number counts (|Rozo et al.ll2007h . 

The primary richness estimator used here is A200, de- 
fined as the number of E/S0 ridgeline cluster members 
brighter than 0.4 L* (in ^-band) f ound within r^oo of the 
cluster center (Han sen et al.l [20051 ) . Analysis of the mean 
weak lensing profiles of th e clusters stacked in bins of 
A200 (|Johnston et al.l 120071 . and § 12. 3[1 provides an im- 
proved measure of the scaled apertures, T2oo- These radii 
are smaller, by ~ 50% on average, than r^o- Our mass and 
luminosity estimates in this letter consistently employ V200 
derived from the weak lensing profiles. However, we remain 
consistent with the original catalog's richness parameter, 
A200, defined within the larger r^g radius. 

As a cross-check, we also bin the sample using a sec- 
ondary richness estimator, L200 , the total fc-corrected i-band 
luminosity of the cluster members. The order parameters, 
A200 and L200, are strongly correlated, especially for the 
richest clusters. We use six richness bins with A200 1? 9, 
described in T able U based on those used in the weak lens- 
ing analysis of ISheldon et al] l|2007l ). The total number of 
clusters is 17335. 



2 INPUT DATA 

Input data for this study comes from two large a rea surveys, 
the Sloan Digital Sky Survey (SDSS: lYorkll2000l) and RASS. 
SDSS imaging data are used to select the clusters and to 
measure the weak lensing shear around their centers. RASS 
data provide 0.1-2.4 keV X-ray luminosities, Lx, measured 



1 Here, M200 is the mass within a sphere of radius r2oo that en- 
compasses a mean density of 200p c (z), with p c the critical density. 



2.2 Measuring Lx 

We measure the mean 0.1-2.4 keV X-ray luminosity, Lx, of 
the clusters in each bin by stacking photons from RASS. 
The typical RASS exposure time for maxBCG clusters is ~ 
400 s, too short to allow significant detections for all but the 
brightest individual clusters. By binning on A200 or L200 we 
can take advantage of the large number of maxBCG clusters. 
The stacking method is described in detail in R07, but there 
are a few differences in the analysis for this letter which we 
highlight here. 

We stack all the photons in a scaled aperture from clus- 
ters in a given richness or optical luminosity bin, treating 
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Table 1. (L x \N 20 o) and (M 200 |JV 20 o> 



galaxies, in the absolute shear calibration, and in the cali- 
bration of the mis-centering model. 



./V200 Range 


f200 

{h- 1 Mpc) 


L x [0.1 - 2.4keV] 
(10 42 h~ 2 ergs' 1 ) 


M 200 
(10 14 h- 1 Mq) 


9-11 


0.49 


2.72 ±0.35 


0.441 ± 0.080 


12 - 17 


0.55 


5.65 ±0.37 


0.600 ± 0.085 


18 - 25 


0.64 


12.8 ±0.98 


0.96 ±0.13 


26 - 40 


0.77 


30.6 ±2.3 


1.68 ±0.23 


41 - 70 


0.90 


56.7 ±4.5 


2.52 ±0.36 


71 - 188 


1.20 


209 ± 31 


5.69 ±0.88 



the brightest cluster galaxy selected by the maxBCG algo- 
rithm as the center of each cluster. In R07 we show that 
this assumption for centering does not introduce a signif- 
icant bias in the calculation of Lx- For consis t ency with 
the cluster mass calculations in Ijohnston" et alJ (|2007ft . we 
scale all radii and luminosities to z = 0.25. This scale red- 
shift is slightly larger than used in R07, but does not sig- 
nificantly change Lx- After building scaled photon maps, 
we construct images, radial profiles, and spectra, and cal- 
culate an unabsorbed, rest- frame mean luminosity, Lx, for 
each bin. The X-ray luminosities are measured within scaled 
apertures, r2oo- These radii are smaller than those used in 
the analysis of R07. However, as Lx is proportional to the 
gas density squared, most of the flux originates near the 
dense core and the effect of this re-scaling is modest, reduc- 
ing the measured values of Lx by ~ 10% — 15%. 

The best-fit spectral temperatures are used to fc-correct 
the observed 0.1-2.4 keV emission to the rest frame at z = 

0. 25. Values for the case of A200 binning are listed in Table 

1. As shown in R07, the data are well fit by a power law 
form, 



{L X \N 200 ) = (42.1 ± 1.7) (JVaoo/40) 1 



.. , n 42 -1 

x 10 erg s 

(1) 



2.3 Measuring M 20 o 

The mean projected surface density contrast, measured from 
the shear of faint background galaxies as a function of pro- 
jected distanc e from the cluster center was calculated by 
ISheldon et al.l (|2007l ) for the binned samples employed in 
this letter. However, as maxBCG clusters are not guaran- 
teed to be centered on the mass density peaks of dark halos, 
Ijohnston et al.l (|2007l ) calculate mean cluster masses using a 
mass model that incorporates a degree of halo mis-centering. 

Specifically, the projected surface density contrast in a 
given bin is modeled as a sum of four components: i) a small 
scale component from the stellar mas s of the central galaxy; 
ii) a mean NFW (|Navarro et al.|[l997f ) halo; iii) a large-scale 
component arising from the clustering of halos (the two-halo 
term); and iv) a correction for the subset of clusters with 
BCG centers offset from halo centers. The mean halo mass, 
M200 , and its uncertainty are obtained by marginalizing over 
the posterior probability distribution using a Markov Chain 
Monte Carlo (MCMC) method. The resultant masses, listed 
in Table [TJ scale as a power-law in richness, 

(M 20 o|iV 2 oo) = (2.1±0.3) (A' 2 oo/4O) 1 ' 28±o ' O4 xiO 14 Af , (2) 

where the error in the intercept is dominated by system- 
atic uncertainties in the photometric redshifts of background 



3 THE Lx-Ahoo RELATION 

In Fig. Q] we plot the mean luminosities and lensing masses 
using binning in -/V200 (Table [1] filled symbols) or in L200 
(open symbols). The two binning choices offer consistent 
results. The gray band shows the best-fit maxBCG-RASS 
relation (±ler) from N200 binning, 



(Lx|iV2oo> = 12.6±i;| 



(MgoojiVgoo) 
10 14 h- 1 Mq 



1.65±0.13 



1 n42 

10 ergs 



(3) 

We estimate an additional systematic error of ±1.6 in nor- 
malization due to the mass scale uncertainties described in 
the previous section. The two additional lines show results 
from X-ray flux-limited cluster samples. The dashed line 
shows the results from the HIFLUGCS clusters of RB02 us- 
ing hydrostatic mass estimates. The data has been re-fit by 
SOB after converting their luminosity and mass estimates to 
ACDM cosmology (see § 3.4 in S06). The dot-dashed line is 
the X-ray count-matching result of S06 for a ACDM model 
with Q m = 0.24, spectral index n 3 = 1, and power spec- 
trum normalization as = 0.85; the inferred scatter in L x 
at fixed mass is (T\ u l\m = 0.40. Both the RB02 and S06 
results are quoted at the median redshift of the maxBCG 
cluster sample, 2 = 0.25, by assuming self-similar evolution, 
L ~ pi {z). This raises the Lx normalization of each by 
27%. Finally, to properly compare with the algebraic mean 
luminosity of maxBCG-RASS, we adjust the normalizations 
of RB02 and S06 by a factor exp(0.5of ni | M ), resulting in an 
~ 8% increase in luminosity at fixed mass. 

The slope of the maxBCG-RASS relation, 1.65 ± 0.13, 
agrees well with the values found previously. The RB02 re- 
sult, 1.50±0.08, is somewhat shallower, but not significantly 
so, and the S06 result, 1.59 ± 0.05, is consistent with the 
present determination. The normalization of the maxBCG- 
RASS relation is consistent with the S06 results, with the 
RB02 masses lying 2a low. Further discussion of these re- 
sults is in §[4] below. 

The maxBCG-RASS normalization can be compared 
with recent lensing results pub lished for smal l samp les of 
massive cluster s at z ~ 0.25. iBardeau et al.l (|2007h and 
iHoeks tra (2007) performed weak lensing analyses of 11 and 
20, respectively, individual clusters. Both studies focus on se- 
lected high luminosity clusters, with Lx > 10 44 /i -2 ergs -1 , 
at redshifts near the median of the maxBCG catalog. Both 
papers find a correlation between Lx and M200, but with 
noisy lensing mass estimates. As each study employs a 
different convention for Lx, we have recalculated Lx for 
these clusters using RASS photon data (see § 4.1 in R07). 
Fixing the slope at 1.6, we compute the normalization of 
L x /10 44 fe~ 2 erg s" 1 at 10 15 fe" 1 M n fin ding 11 ± 2 and 
5.0 ± 0.6 for the lBardeau et all (120071 ) and lHoekstral (120071) 
cluste rs, respectively. At this high mass scale, the Hoekstral 
(120071) normali z ation is consistent with our value, while the 
IBardeau et al l (|2007l ) is ~ 2a larger. A full accounting for 
this discrepancy is beyond the scope of this work. 
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Figure 1. Points show maxBCG-RASS (algebraic) mean Lx and M200 values found by binning on -/V200 (solid circles) or L200 (empty 
circles). The dark gray band represents the ±lcr contours on the best fit relation using the A^oo bins. The dot-dashed line shows the 
S06 relation, while the dashed line shows the S06 fit to the HIFLUGCS clusters from RB02 based on hydrostatic masses. Both relations 
have been scaled assuming self-similar evolution (see §|3J. The error bar in the legend shows the typical la systematic error in the SDSS 
lensing masses, representing an overall shift in normalization that is possible in the maxBCG-RASS relation. The inset plot indicates 
the effects on the slope due to covariance between Lx and -/V200 a t fixed mass (r), as described in §|4] If Lx and A^oo at fixed mass are 
correlated (r = 0.7), this will bias the slope steeper, and if they are anti-correlated, this will bias the slope shallower. 



4 DISCUSSION 

The good level of agreement displayed in Fig.[T]among three 
independent approaches to determining the Lx~M relation 
indicates that optical and X-ray selection methods are find- 
ing similar populations of massive halos. Furthermore, the 
maxBCG-RASS result extends to a lower mass scale than 
is probed by RB02 and S06. In this section, we point out 
effects that could lead to differences among the three mea- 
surements. The discussion is aimed at raising issues to be 
addressed by more detailed analysis in future work. 

Non-zero bias in hydrostatic mass est imates, displa yed 
in early, low-resolution gas simulations l|Evrardl ll99Cll ). is 
a possible source of systematic error that would shift the 
RB02 result relative to the true relation. Recent studies us- 
ing mock X-ray exposures of numerical simulations predict 
a systematic un derestimate of binding mass at the l evel of 
-0.25 in InM jRasia et al.l 120051 ; iNagai et al.ll2007h . Cor- 
recting the RB02 result by this amount would more closely 
align it with the maxBCG-RASS relation. Assuming the 
latter is an unbiased estimate of the underlying halo rela- 
tion, the luminosity offset between the two relations mea- 
sures the Malmquist bias arising from the X-ray flux limit 
of the HIFLUGCS sample used by RB02. Good agreement 
would signal a small bias, meaning small intrinsic scatter 
(< 10%) between luminosity and mass. Such small scatter 
is considered unlikely by the analysis of S06. 

A separate argument can be made based on slope esti- 
mates. If hydrostatic mass estimates scale with true mass as 
(M cst ) oc M t ^l, then one would expect the RB02 slope to 
differ by 1.5e from the maxBCG-RASS value. The measured 
slope difference, 0.15 ± 0.15, implies e = 0.1 ± 0.1. Strongly 
mass dependent hydrostatic biases are therefore ruled out. 

One could shift the RB02 result to higher masses with- 
out requiring a major reduction in scatter, a\ n M\L, as con- 
strained in S06. This would require shifting the S06 result 



by modifying the assumed cosmology. The luminosity nor- 
malization is sensitive to power spectrum normalization, 
Lx ~ °"g^ 4 , so raising erg to 0.95 would shift the S06 result to 
lower Lx and preserve the current level of Malmquist-bias 
for the RB02 result. However, this adjustment would offset 
the SOB and maxBCG-RASS relations at the 2a level. 

Since it is binned by richness, the maxBCG-RASS 
result is sensitive to covariance among Lx and N200 at 
fixed M200. Simulations suggest mild anti-correlation, as 
at fixed mass high concentration halos have higher Lx 
but fewer galaxies (Wechslc r et alj 120061 ). As an illustra- 
tion of the effect of covariance, consider the case of a bi- 
variate, log-normal distribution for Lx and A^oo with con- 
stant covariance. The off-diagonal term can be character- 
ized by the correlation coefficient, r = (<5inz,<5injv), where 
8\nx = QilX— lnX)/a\ n x are the normalized deviations from 
the mean relation. 

Consider a mass function that is a local power-law, 
dn/dlnM ~ M~ a = e~ QM , where ^ = InM. Convolving 
this function with the bivariate log-normal, and using Bayes' 
theorem, allows one to write the conditional likelihood 
P(£,fi\v), where £ = InLx and v = lnA^oo- The result is a 
bivariate Gaussian with mean mass jX{v) = fioiy) — aa 2 ^ u , 
with fio{v) the inverse of the input mean richness-mass re- 
lation and <t m |„ the scatter in mass at fixed richness. The 
X-ray luminosity at fixed optical richness is distributed in a 
log-normal manner with mean 

I(u) = p ■ (p,(y) + ara^a^t) , (4) 

and variance 

a\v = P 2 ■ + <?l\e ~ Zra^cr^e), (5) 

where p is the slope of the halo Lx — M200 relation. 

When Lx and A r 2oo are independent (r = 0), the mean 
luminosity reflects that of the mean mass selected by the 
richness cut. When r ^ 0, the mean is shifted by an amount 
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that scales linearly with the mass function slope a. If the 
scatter is constant, then this shift affects all A^oo-binned 
points equally, and the slope of the richness-binned LX-M200 
relation will be unbiased, d£(y)/dji(y) = p. However, in cold 
dark matter models, the slope a runs with mass, and this 
running can induce a change in slope at the level Ap/p = 

The inset of Fig. [1] explicitly demonstrates this effect 
for the case p — 1.6, <7 M u = 0.25 and — 0.5. We use 
a discrete set of ha los from the Hubble Volume simulation 
|Evrard et al| 2002) that define a mass function which is not 
a power-law Uenkins et al.ll200ll ). To each halo, we assign 
richness and luminosity using a constant log-normal covari- 
ance, then bin the sample on optical richness and compute 
mean luminosities and masses for each bin. Imposing a cor- 
relation coefficient r — ±0.7 tilts the richness-binned Lx~ 
M200 slope by ~ 0.2, or roughly 1.5cr given the error in 
Eq. ©. 

As discussed in iNord et all <|2007h . covariance can both 
tilt the scaling relation and modify the scatter. The second 
observable quantity affected by r is the variance in Lx at 
fixed richness N200, Eq. ©. This scatter was constrained 
in R07 to be oi,\ v = 0.86 ± 0.03. In our model, it takes 
on values of 1.1, 0.90, and 0.72 for r = -0.7, 0, and 0.7, 
respectively. We stress that these values are sensitive to the 
assumed degree of mass scatter at fixed A^oo- Our assumed 
value of 0.5 is smaller than that inferred by the spread in 
galax y velocity dispers ions at fixed ri c hness dBecker et al.l 
120071 ) and employed by I Johnston et all (|2007f) in the lensing 
analysis. Using this larger scatter would result in o t \ v > 1 
for all r 6 [—1,1]. If the R07 measurement of ot\ v < 0.9 
is correct, then either the lBecker et al.l (|2007T ) scatter is an 
overestimate or the simple log-normal model employed here 
is insufficient. 

Improved understanding of the variance is clearly de- 
sired. Constraints on oi\ v can be obtained by individual 
X-ray follow-up observations of richness-selected maxBCG 
subsamples. Such data can also generate hydrostatic mass 
estimates, enabling constraints on the mass scatter er M i„. 
Such a study would begin to test all the elements of the full 
covariance P(£, v\[i) used to describe the underlying massive 
halo population. 



5 SUMMARY 

In this letter, we present a new measurement of the scal- 
ing between X-ray luminosity and mass using ~ 17, 000 
optically-selected clusters of galaxies. X-ray luminosities are 
derived from RASS photon maps, while masses are deter- 
mined independently through modeling of the mean cluster 
lensing profiles. Binning clusters by optical richness, we find 
that the mean Lx and M200 values follow a power-law re- 
lation with slope 1.65 ± 0.13. The consistency in slope with 
previous studies, including those based on hydrostatic mass 
estimates, confirms that the optically-selected maxBCG cat- 
alog selects a population of massive halos similar to those 
of X-ray samples. The consistency implies that systematic 
errors in hydrostatic mass are not strongly scale-dependent. 
Furthermore, we demonstrate how the slope can be affected 
by covariance in Lx and A^oo at fixed mass, and show that 
this covariance is not large. 



Follow-up X-ray observations of -/V200 selected subsam- 
ples would provide key information on the covariance, as 
expressed in Eq. (J5J). Data from future deep surveys, se- 
lected via optical, X-ray and Sunyaev-Zel'dovich signatures, 
will provide a much richer environment for addressing the 
complex interplay between cosmology and covariant scaling 
relations. 
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